Symmetry and transport property of spin current induced spin-Hall effect 
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We study the spin current induced spin-Hall effect that a longitudinal spin dependent chemi- 
cal potential qVs^x,y.z induces a transverse spin conductances G"" . A four terminal system with 
Rashba and Dresselhaus spin-orbit interaction (SOI) in the scattering region is considered. By using 
Landauer-Bitttiker formula with the aid of the Green function, various spin current induced spin- 
Hall conductances G"" are calculated. With the charge chemical potential qVc or spin chemical 
potential qVs=x,y,z, there are 16 elements for the transverse conductances Gp" = Jp,f_i/Vv where 
= x,y,z,c. Due to the symmetry of our system these elements are not independent. For the 
system with G2 symmetry half of elements are zero, when the center region only exists the Rashba 
SOI or Dresselhaus SOI. The numerical results show that of all the conductance elements, the spin 
current induced spin-Hall conductances G"" are usually much greater (about one or two orders of 
magnitude) than the spin Hall conductances G"^ and the reciprocal spin Hall conductances . So 
the spin current induced spin-Hall effect is dominating in the present device. 



Introduction: Recently, an interesting phenomena, the 
spin Hall effect, has been discovered in the spin-orbit in- 
teraction (SOI) system and has attracted considerable 
attention. In this effect, a longitudinal external bias 
(or named charge bias hereafter) or electric field induces 
a transverse spin current or spin accumulations along 
transverse edges. The spin-Hall effect can either be ex- 
trinsic or intrinsic. The extrinsic spin-Hall effect is due 
to the spin dependent scattering[l| and has been found a 
few decades ago. On the other hand, the intrinsic spin- 
Hall effect is due to the SOI, predicted by Murakami 
et.al. and Sinova et.al. in a Luttinger SOI 3D p-doped 
semiconductor and a Rashba SOI 2D electron gas,[l| 
respectively. Since then, a great deal of sequent works 
have focused on this interesting effect. On experimen- 
tal side, two groups by Kato et al[l] and Wunderlich et 
al. have observed the transverse opposite spin accumu- 
lations near two edges of their devices when the longitu- 
dinal voltage bias is added. More recently, a third group 
by Valenzuela and Tinkham took the electronic measure- 
ment of the spin Hall effect, f6l| and they have observed an 
induced transverse voltage in a diffusive metallic conduc- 
tor when a longitudinal net spin current flows through it. 

Very recently, the reciprocal spin-Hall effect [tI, HI has 
been investigated, in which a transverse charge conduc- 
tance is induced by the driving of a longitudinal spin de- 
pendent bias (named spin bias hereafter). The Onsager 
reciprocal relationfzil between the spin-Hall effect and its 
reciprocal phenomenon has been found, and the spin- 
Hall conductance is predicted to be equal to reciprocal 
one. So far, the spin polarization direction considered in 
most of papers on the spin- Hall effect [5, and its recip- 
rocal effect, is along z direction that is perpendicular 
to the 2D electron gas. Since the spin is a vector and it 
can also be polarized in the x, y, as well any n directions. 
So the direction of the output spin in the two transverse 
terminals (i.e. the terminal 2 and 4 in Fig.l) can be in 
the plane of the 2D electron gas, e.g. the x or y direc- 



tion, which has been studied by Nikolic et al. [1, [13] The 
results show that while under a longitudinal charge bias, 
the spin currents of the y direction in the two transverse 
terminals are equal, i.e. they are simultaneously flowing 
out or flowing in (see the left inset in Fig.l). On the 
other hand, in the reciprocal spin-Hall effect, the spin 
polarization direction for the spin bias in the two lon- 
gitudinal terminals (i.e. the terminal 1 and 3 in Fig.l) 
can also be in the plane of the 2D electron gas. A few 
theoretical works have studied this case, and found 
that the charge currents in the terminal 2 and 4 simulta- 
neously flow out or in (see the left inset in Fig.l), while 
the spin polarization of the spin bias is in the y direction. 
This result is different from the usual reciprocal spin-Hall 
effect with the z-direction spin bias, in which the charge 
current flows in at one transverse terminal and out at the 
other one (see the right inset in the Fig.l). 

To explore the vector nature of spin-Hall effect, we set 
the longitudinal bias to be spin-dependent and examine 
the transverse spin conductance, namely the transverse 
spin currents are induced by the longitudinal spin bias, 
which is named the spin current induced spin-Hall effect 
hereafter. So far most of studies focus on the spin-Hall ef- 
fect and its reciprocal effect, less attention has been paid 
to the spin current induced spin-Hall effect. It is the pur- 
pose of this paper to fill up this gap. We study the spin 
current induced spin-Hall effect in a finite mesoscopic sys- 
tem by using Landauer-Biittiker formula with the aid of 
the Green function. Considering a four-terminal ballistic 
2D rectangular region and the center rectangular region 
having the Rashba and Dresselhaus SOI (as shown in 
Fig.l). A spin bias Vs (s = x, y, and z) is added at 
the longitudinal terminals 1 and 3.jT3| The spin bias can 
be provided from the device of the spin cell, which has 
been suggested by some recent theoretical works lla and 
also been realized in a few experimental works. [lj| Un- 
der the spin bias K, the chemical potentials for the 
spin-up and spin-down electron in the terminals 1 and 3 
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FIG. 1: (color online) Schematic diagram for the four termi- 
nal rectangular sample with the Rashba and Dresselhaus SOI 
in the center region. The four leads are ideal without SOI. 
The spin bias Vs or the charge bias Vc are added on the lead- 
1 and lead-3, the induced transverse spin current Jp,s or the 
charge current Jp,c are probed in the lead-2 and lead-4. The 
left and right insets depict the directions of the transverse 
spin or charge currents when the currents are simultaneously 
flowing out (G2 = G4) and the currents are flowing in at one 
terminal and out at the other one (G2 = — G4), respectively. 



are split, and /zi,sT = -Mi.si = -Ms.sT = M3,si = 6^/2, 
here s = x, y, and z represents the spin polarization 
direction. For comparison, we also consider the charge 
bias Vc on the longitudinal terminals 1 and 3, in this 
case = fj,i^si = -M3,sT = -f^s.si = The 

transverse two terminals 2 and 4 act as the measuring 
terminals, and their chemical potentials (or the termi- 
nal bias) are set to zero. Due to the SOI, the spin bias 
Vs or the charge bias Vc can induced a transverse spin- 
Hall currents Jg as well as the charge-Hall current 
which will be investigated in this paper. For this pur- 
pose, we use four component vectors J = {Jx, Jy, Jz, Jc) 
and V = {Vx,Vy,Vz,Vc) to represent the spin (charge)- 
Hall current and spin (charge) bias, respectively. In the 
small bias limit, the relationship between J and V can be 
characterized by a 4 x 4 matrix C^" (p, v G {x, y, z, c}). 

We have studied the relations among these matrix el- 
ements from the symmetry point of view. Four cases 
have been discussed: (1). only Rashba SOI is present, 
(2). only Dresselhaus SOI is present, (3). both of the 
Rashba and Dresselhaus SOI are present, (4). the Rashba 
and Dresselhaus SOI have the same interaction strength. 
The results are summarized in the Table I. All elements 
are found to have the property \G*i,^\ = |G4''|, where 
G*^" = Jp.p/K with j3 = 2,4.[ll| This means that the 
absolute value |J2.^i| of the spin (/x = x,y,z) current or 
the charge = c) current in the terminal 2 is equal 
to one 1^4,^1 in the terminal 4 regardless of longitudinal 
driving bias V^. From Table I we see that if only the 
Rashba SOI or the Dresselhaus SOI is present in the cen- 
ter region, half of the matrix elements of the transverse 
conductances C*"^ are zero. We have also numerically 
studied the behavior of these non-zero spin matrix ele- 
ments and find that the spin current induced spin-Hall 
conductances are much larger (about one or two orders 



of magnitude) than the spin Hall conductance as well as 
its reciprocal one. 

Hamiltonian and Solution: The system that we con- 
sidered is sketched in the Fig.l. The central gray region 
is the semiconductor sample in which the SOI of Rashba 
or/and Dresselhaus are present. The sample is connected 
to four ideal non-SOI leads. The Hamiltonian of the cen- 
tral region is given by Hq — /2m* + V{x,y) + j- {a^Py ~ 

o'yPx) — ^{o'xPx — o'yPy), whcrc a and f3 are the coefficients 
of the Rashba and Dresselhaus SOI, V{x, y) is the hard 
wall confining potential. In the tight-binding representa- 
tion, the total Hamiltonian can be written as: [13, HI, [3| 

H = y^\-ta\ai+sx - ta\a.i+sy + ia\{VR(7y + VD(Jx)ai+Sx 

i 

- ia\{VR(Tx + VDiyy)ai+Sy] + H.C (1) 



where a\ — {a[^,a\^), i is the lattice site in the central 
region and the leads, t = h? /2m*a?' is the hopping matrix 
element with the lattice constant a, 5x and 5y are the unit 
vectors along the x and y directions, and Ux , (Jy are Pauli 
matrices. Here Vr — a/2a and Vd — P/2a represent the 
strength of the Rashba and Dresselhaus SOI, respectively. 
Vr and Vd are non-zero only in the central gray region. 

Since there is no SOI in the leads, the particle cur- 
rent Jp^sa in the lead-p (p=2,4) with spin at ±s direc- 
tion (cr ~] or I stands for the -\-s or — s direction with 
s = x,y, z) due to the longitudinal spin or charge bias Vg, 
can be obtained from the Landauer-Biittiker formula: [20j 

Jp,sa — (-^/^) ^^(7.(t' '^psa,qs' <t' {^p.sa l^q,s'cr') whcrC fip^sa 

and Hq^s'a' are tne spin-dependent or spin independent 
chemical potential related to the bias VJj , which has been 
detailed in the the introduction. Tqsia',psiy is the trans- 
mission coefficient from the lead-p with spin sa to the 
lead-g with spin s'a' . Note that the spin index f and J, 
in the Hamiltonian (1) represent the spin pointing to the 
+z and — z-direction, and not pointing to the ±x and 
±y-direction. So in order to calculate the transmission 
coefficient Tqs'a-',ps<7 for s 7^ z and/or s' ^ z, we need 
to rotate the z-axis in the spin space in the lead-p (lead- 
q) to the s-direction (s'-direction) by taking an unitary 
transformation: 



s(s') 



flit 
an 



(2) 



where i is the lattice site in the \ead-p{q). The operator 
flio- in the center region and other two leads does not 
change. The unitary matrix in the Eq.(2) is 



cos 7 



(3) 



where 



is the directional angle of s-direction and 



U^CTsUs — Uz- Under this unitary transformation, the 
Hamiltonian of the leads does not vary because that the 
leads' Hamiltonian does not include the Pauli matrix, 
only the Hamiltonian that describes the coupling between 
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the leads and the center region changes. After the unitary 
transformation, the z-axis of the spin in the lead-p(g) is 
in the s(s')-direction, then the transmission coefhcient 
Tqs'a'.psa Can be easily obtained as:[20j 

T^qs' <j' ,ps(7 — Tt\Y^ qg' fj'G s^G ], (4) 

where the line-width function Tp,scr = — = 

Tp ® UsScrUj, the retarded self-energy Sp^^ = Sp 
UgSo-Ul with Tp and Sp, respectively, the line-width 
function and the self-energy of the lead-p for no-spin- 
index lattice system. Here the matrix is: 



St = 













2> 







(5) 



The retarded and advanced Green's functions C'" in 
the Eq.(4) can be calculated from ^ G'' = [G"]t = 
{EpI—Hq—^ ^ Sp^^}~^, with the unit matrix I. Since 



X;^(U.S<.Uj) = ® I that is inde- 



pendent of spin-direction index s, the Green functions 
C"'" are also independent of spin-direction index. This 
in turn indicates that the Green functions C'" remain 
unchanged in the above unitary transformation, in which 
the spin axis in the leads is rotated. 

Once the particle currents Jp^sa [p = 2, 4) are obtained, 
the spin current Jp,s and the charge current Jp.c can be 
obtained straightforwardly: Jp,c — ^{Jp,s\ + Jp,si\ and 
Jp,s = (^/2){Jp,sT ~ Jp:si\ where s = x,y,z. In fact, it is 
easy to show that the charge current Jp c is independent 
of spin-direction index s of the lead-j> (i.e. Jp,x\ ~^ 'Jp,xi — 

^p^yl ~^ ^p^yl — '-^p^^\ ~^ '^P:^l — ^p^^^l ~^ ^p,^l^: the 
subscript s is neglected from now on. Then from the 
current Jp ~ iJp,x, Jp,y, Jp,z, Jp,c) and its driving bias 
V = {Vx,Vy,Vz,Vc), we get all the sixteen conductance 
matrix elements G^"^ = Jp^^i/Vi, {fi,i^ G c. x. v. z).[l6j 
They are the charge Hall conductance Gp'^ describing the 
transverse charge current Jp^c induced by the longitu- 
dinal charge bias Vc, the spin Hall conductance Gp'^ for 
the transverse spin current Jp^s induced by the longitudi- 
nal charge bias Vc, the reciprocal spin Hall conductance 
Gp^ representing the transverse charge current Jp.c in- 
duced by the longitudinal spin bias Vg , and the spin Hall 
conductance Gp* induced by spin current describing the 
transverse spin current Jp^g induced by the longitudinal 
spin bias Vg. 

Relations due to symmetry: We now study the rela- 
tions among the sixteen conductance matrix elements 
G^" considering the symmetry of the device. First, our 
system, satisfies the time reversal symmetry, [2l| so the 
transmission coefficients satisfy the relation: Tpg„^qgi„i = 



j., which determines the properties of the trans- 



verse conductances. Second, if the shape of the device 
has the geometrical symmetry, the results can be greatly 
simplified and many matrix elements are zero. In the 
following, we consider the rectangular center region that 
has G2 symmetry. 

For the rectangular sample, the shape of the device [i.e 
the confining potential V{x, y)] is invariant under the ro- 
tation transformation G2 by rotating an angle tt around 



the a;,?/, z-axis at the center point. But the SOI part 
of the Hamiltonian is varied under the space rotation 
transformation G2 . In order to keep the invariance of the 
total Hamiltonian iJ, we construct the unitary transfor- 
mation U = Uo (S^ Us, where Uq and Ug are the rotation 
transformation in the real space and in the spin space, 
respectively. In the following, we list all of the unitary 
transformations U, under which the total Hamiltonian 
H is invariance (i.e. WHU = H). We then derive the 
relations among the transverse conductance matrix ele- 
ments from the symmetry. We consider four cases of SOI: 
(1). with Rashba SOI only, (2). with Dresselhaus SOI 
only, (3). with both the Rashba and Dresselhaus SOI, 
(4). similar to (3) but the strengths of the Rashba and 
Dresselhaus SOI are equal. To illustrate the derivation, 
we shall discuss the case of the Rashba SOI case in detail. 
For other three cases, similar discussions apply and we 
only give the results. 

First we study the case with only Rashba SOI using 
the following transformations. 

(i) Ui = C2Z s.xp(—i^a z) , which performs the trans- 
formation: X —>■ —X, y — > —y, z ^ z, and ax — > —(^x, 
Uy — > —<Jy, <Jz — > fz- Under this transformation Ui, the 
system (including both of the real space and the spin 
space) is rotated by an angle tt around the z-axis, which 
leads J2,xa J4,xa and fii^xa ^^3,xa, so the spin cur- 
rent J2^x = {fl/^){J2,x^-J2,xl} {h/2){J4^xl-J4,x^} = 

~J4:,x, and the spin chemical potential Vx remains un- 
changed because the charge chemical potential Vc and 
the <7x change signs due to the rotations in real space 
and spin space, respectively. Similarly, the transforma- 
tion for the others J2.fj, and can be also deduced as: 



(6) 



Then the relations among the matrix elements of the 
transverse conductance can be obtained straightfor- 
wardly. For example, Gf^ — J2.x/Vx = —Ja.xIVx = 
-Gf^, etc. 

(ii) C/2 = C2y®ex-p{—i'^ax), which performs the trans- 
formation: X — > ~x, y ^ y, z ^ —z, and <Jx — > <Jx, 
ay —o'y, az — > — CTz. This transformation U2 is equiv- 
alent to rotate the real space by tt around the y-axis and 
to rotate the spin by tt around the a;-axis. Under the 



J2.,c - 


Ji,c 


Vc - 


-Vc 


J2,x - 


~Ji,x 


Vx ~ 


Vx 


J-2,y - 


—Ji,y 




- Vy 


J2,z - 




Vz - 


^ -Vz 



unitary transformation U2, 
transformed into: 



'Jp,!^ (P = 2,4) and V^ are 



J. 



2(4)x 
h{4),x 
^2(4), y 
h(4),z 



J. 



2(4)x 
h{4),x 



-J'. 



-J'. 



2(4),a 



2(4),z 



Vc 
Vx 



'Vc 
-Vx 



Vz 



(7) 



Then Gf = Jp.c/K = Jp,c/(-K) = -G^^ 



so Gf = 



(p = 2,4). In fact, from u\HU2 = H, we obtain that 
eight matrix elements of the transverse conductance are 



zero, G'^ = 
[p - 



G 



Gcx 
p 



2,4). 



Gxc 
p 



Gxx 
p 



GVpy 



Qyz 
\-j p 



G^y 
\-j p 
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TABLE I: Symmetry of the transverse spin or charge con- 
ductance for the system with the Rashba and/or Dresselhaus 
SOI: (a), Vr / 0, = 0; (b), Vr = Vd ^ 0; (c), 
Vr^Vd ^ Q and (d), Vr = Vd^ 0. The symbol '0' in- 
dicates the corresponding G^"^-^ = 0, and the symbol '+ (— )' 
denotes Gj'^ = ±64". 
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(iii) U3 = C2x'Siexp{—ijay), which performs the trans- 
formation: a; — s- a;, y — > —y, z —z, and — > — cr^, 
ay —^ ay, az — > —az- This transformation U3 rotates 
the real space by tt around the x-axis and the spin by 
TT around the y-axis. Under this transformation, Jp_^ 
{p = 2,4) and are transformed into: 



J2,c - 


Ji,c 


Vc - 


V, 


J2,x - 


~* ~Ji,x 


Vx - 




J2,y - 


Ji,y 


Vy - 


- Vy 


Jl.z - 


—Ji.z 


Vz - 


-Vz 



(8) 



It is worth mentioning that the above three unitary trans- 
formations are not independent on each other, and we 
have U3U1 = U2- So Eq.([7|) can be obtained from Eq.® 
and Eq.®. 

Combining Eqs. (|6l7l8p . the relations among the ma- 
trix elements of the transverse conductances G^'^ are ob- 
tained and summarized in the Table 1(a). The sixteen 
conductance matrix elements are arranged in the ma- 
trix form, in which the columns (labelled by J^) denote 
the measured transverse conductances (from lead 2 or 
4) and the rows are the longitudinal bias Vi, in lead 1 
or 3. The zero conductance matrix elements are indi- 
cated by the symbol '0'. The symbol '-I- (— )' denotes the 
non-zero conductance matrix elements G^", which have 
the same (opposite) signs for the lead-2 and lead-4, i.e. 
Gj" = ±Gl^ . For some pair of conductance matrix ele- 
ments, they may have the same value or differ by a sign, 
e.g. Gp^ = Gp"^, p = 2,4, which is marked by the let- 
ter symbols in the Table I. Of which, different symbols 
present different conductance values, and the symbol a 
(or a) are for G^^ = a (or G2" = —a). 

To examine the system with Dresselhaus SOI only, 
we need only change ax to ay and ay to ax in the 
above Hamiltonian with Rashba SOI and unitary trans- 



formations. As a result, the above three transforma- 
tions are changed to Ui — C^z ® exp{—i^az), C/2 = 
C2y <8) exp{—i^ay) and U3 = C2X €5 exp{—i^ax)- Under 
these transformations, we find all the symmetry relations 
among the sixteen conductances and the results are listed 
in the Table 1(b). For the case of having both the Rashba 
and Dresselhaus SOI {Vr ^ and Vd 7^ 0), only the 
transformation Ui can keep the Hamiltonian invariant, 
and it leads the symmetry relations shown in the Table 
1(c). For the case of Vr = Vd ^ 0, except Ui, there is an 
additional symmetry U4: I ® (cx^y) cfy{x)i'^z ^ (^z), 
which leads to the symmetry between the Sx and Sy, and 
the relations of the conductance elements G^'^ are ex- 
pressed in the Table 1(d). In fact, Ui only takes the 
transformation in the spin space, so the symmetry Ui 
always exists and G^'^ = G^^ = for any confining po- 
tential V{x,y) (not requiring the rectangular sample). 

From these symmetry relations, we conclude: (i) All 
the conductance matrix elements are found to obey the 
relation \G2"\ = IGjj"'!, i.e., the absolute value of |J2,/j| 
of the spin (/i = x, y, z) or the charge (/i = c) current 
in the terminal 2 is equal to the one of \ Ji,^i\ in the ter- 
minal 4. Furthermore, G^ — —G^ for the eight block- 
diagonal elements in Table I, and G^ = G4'' for the 
eight non-block-diagonal elements, (ii) Exchanging ax 
and ay, the system with Rashba SOI is same as that 
with Dresselhaus SOI, which can be seen from the Table 
1(a) and 1(b). (iii) If only the Rashba SOI or Dressel- 
haus SOI is present in the center region, half of matrix 
elements G^"^ are zero.fl^l (vi) For the usual spin Hall 
effect or its reciprocal effect (in which the spin is polar- 
ized along the z-direction) , the transverse spin current or 

charge current are conserved (i.e. Gj'^*-'^^-' — —G^'^^'^^''). 
Furthermore, G^"^ and G'^^ satisfy the Onsager relation: 
^2(4) = G'2(4)-0 111 While for G'=y and G?'^ as well 
the non-block-diagonal elements in Table I, they either 
flow in or out of the terminals 2 and 4 (i.e. G2 = G4), 
and do not satisfy the Onsager relation, because that the 
Onsager relation requires that the currents in the two 
transverse terminals 2 and 4 must be J2 = — J4. 

Numerical results and Discussion: In the following, we 
numerically study the non-zero conductance matrix ele- 
ments G^'', and we mainly focus the spin Hall conduc- 
tance induced by the spin current G^" (s,s' € x,y,z). 
For simplicity, only the case of Rashba SOI is shown, and 
the results for the other three case have similar physics. 
In the numerical calculation, we consider the square scat- 
tering center region for convenience. But notice that the 
symmetry relation in Table I is valid for the rectangular 
scattering center region. In the calculation the electron 
effective mass m* is set 0.05TOe,[23, HH] and the Fermi 
energy Ep — — 3.8t (in Fig. 4, 5) which is near the band 
bottom — 4i, with t = 1 as a energy unit and the corre- 
sponding lattice constant a « 2.6nm.[l, [13] The size of 
center region L is chosen in the same order as the spin 
precession length Lso {Lso = '^o,t/2VR) over the pre- 
cessing angle tt. If taking Vr = 0.03t (corresponding to 
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FIG. 2: (Color online) The spin Hall conductances C"^ and 
the reciprocal spin Hall conductances C^" vs Fermi energy 
Ef, with the parameters L — 34a and Vr = 0.034. 




Fermi Energy E^ (t) 

FIG. 3: (Color online) The non-zero spin current induced spin 
Hall conductance elements (black), G^^ (green or gray), 
Gl^ (black), and — Gf^ (green or gray) vs Fermi energy Ef. 
The parameters are: L = 34a and Vr = 0.03t. 



a « 1 X 10"" eV-m), Lso « 50a. 

First, Fig. 2 shows the spin Hall conductance G^"^, the 
reciprocal spin Hall conductance C^^, the in-plan spin 
Hall conductance G^^, and its reciprocal conductance 
(^cy they are odd functions of the Fermi en- 

ergy Ef and vanish at the half-filled band with Ef — 0. 
The behaviors of the spin Hall conductance G^'^ and 
G^'^ are similar to the previous works done by Nikolic 
et.al.[^ From the Fig. 2, we can find that out-plane 
matrix elements G^'^ = G"^^, which satisfies the Onsager 
relation. On the other hand, the in-plane matrix ele- 
ments G^'^ and G'^'^ do not obey the Onsager relation, 
because G^'^ and G'^^ are simultaneously flow out of or 



into the two transverse terminals 2 and 4 (as shown in 
the left inset of Fig.l), and they do not satisfy the cor- 
responding condition. G^'^ is smaller than G^*', and both 
of them strongly depend on the Fermi energy. 

Second, we examine the spin Hall conductances in- 
duced by spin current G*" (s, s' £ x,y,z). Fig. 3 shows 
the four non-zero matrix elements G*" versus the Fermi 
energy Ef. Different from G*"^ and G'^^, the spin current 
induced spin Hall conductances G'^^ are even functions 
of Ef, and jG*'* | reaches the largest value at the half- 
filled band with Ef = 0. In particular, the spin current 
induced spin Hall conductances G'^^ are much larger (at 
least one order of magnitude larger) than G*"^ or G"^* (see 
Fig. 2 and Fig. 3). This means that the spin current in- 
duced spin Hall effect is the dominating effect in the SOI 
system. For the spin current induced spin Hall conduc- 
tances G"'* , they do not obey the Onsager relations, in 
general G^*' 7^ G^" and G^^ ^ G^"=. This is because 
that the present device is not simply four-terminal sys- 
tem, more like the eight-terminal system with the index 
of both the terminal and the spin. Also because that for 
the circuit of the Onsager relations, it requires that the 
boundary condition of the transverse terminals 2 and 4 
is J2 — — J4, i.e. in the external circuit, the terminals 
2 and 4, and the terminals 1 and 3, are connected, re- 
spectively. However, in the present system, the current 
^2'^ 7^ — G4'' for the non-block-diagonal elements. On 
the other hand, for the smaU Vr (e.g. Vr < O.OSt), |G^y| 
and |G^^| are still approximatively equal to |G^^| and 
G^^^l, respectively. 

The spin current induced spin Hall conductances G** 
versus the sample size L and the Rashba SOI strength 
Vr are plotted in Fig. 4 and Fig. 5, respectively. With the 
increasing of the size of center region L, the conductances 
G^"^ and G^^ are greatly enhanced (see Fig. 4), because 
that the transverse terminals 2 and 4 are much wider 
at the large L. But at a few special sizes L (e.g. L — 
26,46 for Vr = OMt), G'^y and G^^ are anomalous, in 
which G^y and G^^ are minimum. This is because that a 
subband passes the Fermi energy Ef. The conductances 
G^^ and G^^ have the similar behaviors with G^^ and 
G'-^y while in the small Vr region, i.e. the corresponding 
spin precession length Lso — i^atjlVR, > L. But for 
the large Vr case (i.e. Lso < L), G^^ and G^^ are not 
regular with the size L. Because in the present spin Hall 
device, the quantum states are extended. [2^ 

On the other hand, with increasing of the Rashba SOI 
strength Vr, the conductances G^* are oscillatory (see 
Fig. 5). For the larger Vr, G*** shows a chaotic behavior. 
Furthermore, at the large Vr, |G^^| and |G^^| are not 
equal to |G^^| and |G^^|, i.e. the Onsager relations for 
the spin current induced spin Hall conductances G** are 
completely violated. 

At last, we investigate how the spin current induced 
spin-Hall conductances G"^ are affected by the disor- 
ders. In the previous calculation, the on-site energies 
ei (i.e. in the term X^i^i'^l'^O Eq.(l) are taken to 
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FIG. 4: (Color online) The spin current induced spin Hall 
conductances G"" vs sample size L for different SOI strength 
Vr = 0.03 (a), Vr = 0.06 (b), and Vr = 0.1 (c). The curves 
of Gj^ and G2^, — Gf^ and Gi^ are nearly coincident for the 
small Vr = 0.03 and 0.06. The Fermi energy is: Ef = -3.8t. 
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FIG. 5: (Color online) The spin current induced spin Hall 
conductances G^" vs the Rashba SOI strength Vr, with the 
parameters L — 34a and Ef — —3.St. 



be zero when there is no disorders. To consider the 
effect of disorder, random on-site potentials in the 
center region are added with an uniform distribution 
[—W/2,W/2] with disorder strength W. The conduc- 
tance is obtained by averaging over up to 1000 disorder 
configurations. Fig. 6 shows the normalized conductances 

Gs 



[G: 



SS 



G"'' {W + 0)/C {W = 0)] versus the 



disorder strength W for the different size L. From the 
Fig. 6, we can see that the spin current induced spin Hall 
conductances G"* decrease with the increasing of the 



disorder strength W , but G*^ keep the large value while 
W < t. This behavior is similar to the spin Hall conduc- 
tances which have been investigated recently. [27] Hence, 
in the dirty case, the spin current induced spin Hall effect 
is still dominant in the finite 2D SOI system. In addi- 
tiona, for the different size of the device, the relation of 
^ratio ^ratio vcrsus W almost remain the same. On 
the other hand, Gf^^^^ and G^^jj^ can keep larger value 
in the system with small size than the system with big 
size for fixed disorder strength W. 

Summary: The spin current induced spin Hall effect 
is investigated in a four terminal system with the center 
region having the spin-orbit interaction (SOI). Because 
of the vector nature of spin, the charge and three spin 
components form a 4 x 4 = 16 transverse conductance 
matrix whose matrix elements include the spin current 
induced spin Hall conductances G"^ (s, s' € x, y, z), the 
spin Hall conductances G'^'^, the reciprocal spin Hall con- 
ductances G'^^, and the charge Hall conductance G'^'^. Of 
these matrix elements, we found that in general, G''* are 
much larger (about one or two orders) than the others. 
This means that the spin current induced spin Hall ef- 
fect is the dominating effect in the present device. By 
analyzing the system's symmetry, the relations among 
these conductance matrix elements are found. The re- 
sults indicated that eight matrix elements (^G^^/^v/v^/yy 
and Q^^/^c/c^/cc-^ have conserved quantities as the usual 
Hall effect. But the other eight matrix elements cor- 
respond to the current simultaneously flowing out of or 
into two transverse terminals (as shown in the left inset of 
Fig.l), which is different from the usual Hall effect. When 
only the Rashba SOI or Dresselhaus SOI is present in the 
device, half of matrix elements are found to be zero. 
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